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Abstract. We classify the future of the universe for general cosmological models 
including matter and dark energy. If the equation of state of dark energy is less then 
— 1, the age of the universe becomes finite. We compute the rest of the age of the 
universe for such universe models. The behaviour of the future growth of matter 
density perturbation is also studied. We find that the collapse of spherical overdensity 
region is greatly changed if the equation of state of dark energy is less than —1. 
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1. Introduction 



The Universe is replete with dark energy whose nature is almost completely unknown, 
excepting that its "equation of state", w — p/p, is negative. Dark energy determines 
the future of the Universe. Recently, there is renewed interest in the future of the 
Universe partly because the state of the future of the universe would dramatically be 
changed in the presence of dark energy |I] |2] . Moreover, another motivation comes from 
the fact that the equation of state of dark energy is being constrained by cosmological 
observations and the discovery of a new type of the future singularity with w < — 1 
dark energy (called phantom [HE]): the big rip This singularity is intriguing because 
a spacetime singularity appears even if the weak energy condition (p+p > 0) is violated. 
The big rip is the singularity where the universe expands so rapidly due to the repulsive 
rather than the attractive nature of gravity that both the scale factor and the Hubble 
parameter diverge there. 

In this paper, including w < — 1 dark energy, we classify the future of the universe 
and the conformal diagrams for cosmological models with dark matter and dark energy 
with general equation of state. We also calculate the remaining age of the Universe if 
the age is finite. We also investigate the future evolution of matter density perturbation 
in the universe with dark energy. 

Throughout the paper, we limit ourselves to a constant w. However, the case 
with time varying w could be obtained by combining our results (however, see [7j for 
exceptional cases). 



2. Diagram of the Universe 

The Friedmann equation for the universe with (nonrelativistic) matter and dark energy 
is 




= Q M a- 3 + Q x a- 3 ^ + (1 - Q M - Qx)a~ 2 , (1) 



where Hq is the present Hubble parameter and w is the equation of state of dark energy 
and Qm and fix is present density parameter for matter and dark energy, respectively. 
The scale factor has been normalized such that a = 1 at present. 

The Friedmann equation Eq. Q can be read as the "energy equation" : 

+ V ( a ) = E > ( 2 ) 

where the potential energy term V(a) is 

V(a) = -VlMa- 1 - fljfl- 1 - 3 : (3) 

and the total energy E(= Qk — ~K/Hq) is E = 1 — Qm — Qx- If there are roots of 
E — V(a) =0 with a > 0, the universe can stop its expansion (d = 0) and recollapse or 
bounce. Let us discuss the roots of E — V(a) = in the following. 
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Figure 1. The "phase diagram" of the universe in Q,m — fix plane for w — 
—0.3,— 1/3,— 0.5,— 1,-1.5,— 2. The regions, labeled "no big bang", "expands forever" , 
"big rip", and "recollapses eventually" are divided by the solid lines. The horizontal 
dashed line is Qx = 0. The region of "no big bang" represents the bounce cosmologies. 
In the region of "expands forever" for w > — 1, the universe will expand forever. In 
that of "big rip" for w < — 1, the universe will end in a big rip. The expansion of the 
universe is currently accelerating above the dotted line, while decelerating below it. In 
the region of "recollapses eventually" , the universe will eventually stop its expansion 
and recollapse to a big crunch. 
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A: Qx >0&0>u>> — 1/3. In this case, the dark energy behaves like a matter 
with small negative pressure. If E (= VL K ) < 0, since E — V(a) = 1 > at a = 1 
and E — V(a) — > 0) for a — > oo, there is a root of E — V(a) at a > 1. Hence 
the universe will recollapse in the closed model Qk < 0. Interestingly, if E > there 
is no root and the universe will expand forever. The future of the Universe would be 
quite different from that with a cosmological constant: The Universe with negative dark 
energy density can expand forever without big crunch. 

In top left panel of Fig.l, we show the "phase diagram" of the universe for w = —0.3. 
The two regions, labeled "recollapses eventually" and "expands forever" , are divided by 
solid line. The horizontal dashed line is fix = 0. 



B: w — —1/3. In this special case, the Friedmann equation Eq.((H) is reduced to 
(H/H ) 2 = f2 A /a~ 3 + (l— f2 A /)a~ 2 and is independent ofQx- Dark energy with w = —1/3 
behaves like a curvature, so that the destiny of the universe depends only on the matter 
density. For Qm < 1 the universe will expand forever, while for > lit will recollapse. 
In top right panel of Fig.l, we show the "phase diagram" for w = —1/3. 



C: Qx > O&w < —1/3. The universe can recollapse or bounce if E < V m = 
max(V(a)). Let a m be the scale factor at the maximum of V(a), where a m is given 
by 

-l/3w 

(4) 



n 



M 



3w)tt 



X 



and V m is 



3w Q 



M 



1 + 3w a r , 
Hence the critical condition E - 



Vrr 



is rewritten as, 



1 — Vl m — O 



i-x 



3 w tlx 



-3w 



M 



-1-3™ 



(5) 



(6) 



[1 + 3w) tt x 

This equation is the same as Eq.(5) in Moles (1991) if we set w = —1. The bounce 
would occur if E < V m and a m < 1; the recollapse would occur if E < V m and a m > 1 

In Fig. 1, we display the "phase diagram" of the universe depending on the equation 
of state of dark energy (w = —0.5, —1, —1.5, —2). For w < —1/3, the three regions are 
divided by the two solid lines. The region of "no big bang" represents the bounce 
cosmologies. As w becomes negatively larger, the region of bounce 

becomes larger. This can be understood from 



dw w 



(lna m - ln(fi x /f2 M )) 



(7) 



for fixed Qm and Qx- This implies that for a m < 1 (or Qm < — (1 + 3w)Qx), V m 
increases as w decreases when Om and Qx fixed, so that the universe can bounce more 
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easily. On the other hand, for a m > 1, V m decreases as w decreases when Q M and fix 
fixed, so that the universe can recollapse more hardly. 

If the equation of the state of dark energy is less than — 1 |HEj, then the universe 
expands so rapidly that the scale factor will diverge and the space will eventually 
be torn apart and the universe will consequently result in the "big rip" |B]: the 
singularity with the positively divergent Hubble parameter. In the bottom panels of 
Fig.l (w = —1.5, —2), the region of "big rip" represents the cosmological models with 
big rip singularity. $ 

D: fl x < & w 7^ —1/3. For > w > —1/3 the universe will recollapse if E < V m . 
While for w < —1/3 there is a root at a > 1, since E — V(a) is 1 at a = 1 and 
f2x&~ 1 ~ 3 "' (< 0) at a ^> 1. Hence the universe will recollapse in fix < & w < — 1/3. 

3. The Life of the Universe 

For w < — 1, the age of the Universe becomes finite even if the energy density of dark 
energy is positive. Then the immediate question would be, "how much time is left for 
the Universe?". So we calculate numerically the remaining age of the universe, ti e ft, for 
general non-flat universe with w < — 1: 



where we have introduced x = 1/a. The case of flat universe is calculated in [II 1^ ITU]. 

Fig. El is H t\ e f t for general non-flat universe models: for w = —1.5, H t\ c { t is 1.0 
from top in step of 0.2; for w = —2.0, H ti e{t is 0.6 from top in step of 0.2; for w = —3.0, 
H t\ c f t is 0.3 from top in step of 0.1; for w = —5.0, H t\ e ft is 0.2 from top in step of 0.1. 

Fig. is for a flat model: H t\ e f t is 5.0, 3.0, 2.0, 1.0, 0.5, 0.3, 0.2 from top to bottom. 
From Fig. El we find that the remaining age is primarily determined by fix- Moreover, 
for negatively larger w dependence on Qx becomes weak and the age is essentially 
determined by w (see Fig. 2 for w = —5.0). This fact can easily be found from Eq.(|S|: 
since for \w\ ^> 1 dark energy becomes dominant during most of the remaining age of 
the universe, the integral can be approximately estimated as 



which explicitely shows that ti e f t is more sensitive to w than to fix- It is found from 
numerical calculations that this approximates the integral Eq.fjHJ) well (within 10%) for 



| The properties of the singularity and the existence of a stable fixed point are studied in 






(9) 



0.5<ft x <2.0ifu;<-1.3. 
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Figure 2. ffo^ieft for general universe model: for w = —1.5, -ffoiieft is 3.0 from top in 
step of 0.2; for w — —2.0, 0.6 from top in step of 0.2; for w = —3.0, i?o^ioft is 0.3 from 
top in step of 0.1; for w = —5.0, -ffo^ieft is 0.2 from top in step of 0.1. 



4. Conformal Diagrams 

Similar to the fate of the universe, an interesting question would be "what is the entire 
structure of spacetime for the universe with dark energy?". We show the conformal 
diagrams of flat cosmological models with various dark energy. There are several works 
which have some overlap [12 EH El EH CHI UH UH] • Our aim here is to collect and classify 
these results including matter and w < — 1 case for comparison and for completeness. 
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Figure 3. ff *left for flat models. H ti e ft = 5.0,3.0,2.0,1.0,0.5,0.3,0.2 from top to 
bottom. 

4-1. Dark Energy without Matter 

For simplicity, we first consider cosmological models without matter. The conformal 
diagrams of such cosmological models (for w > —1) are studied in [T2J[T3]. We include 
cosmological models with w < — 1 dark energy as well. 

Since the metric of flat models is conformal to Minkowski spacetime 

ds 2 = -alt 2 + a(t) 2 (dr 2 + r 2 d{l 2 ) = a^f^-dr, 2 + dr 2 + r 2 dtt 2 ), (10) 

the conformal diagram is a subset of that of the Minkowski spacetime and the range of 
the conformal time rj depends on w. 
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A: w > -1/3. Since a(t) = t 2 ^ 1+w \ 

f dt = 3(1 + W) (i +3 ^)/3(i +TO ) 



(11) 



The range of r\ corresponding to that of t (0 < t < oo) is, < n < oo, so that the 
conformal diagram is the upper half of the conformal diagram of Minkowski spacetime. 

B: w = —1/3. Since a = t, n = Int. The range of rj is — oo < f] < oo, and the 
diagram is almost the same as that of the Minkowski spacetime excepting the presence 
of big-bang singularity at the past null infinity. 

Since there is a past null singularity, the physical size of the past light cone dn for 
an observer at r = (particle horizon) is given by 



and it is infinite for w = — 1/3. Hence, there is no horizon problem in this case. 

C: —1 < w < —1/3. Since a(t) = t 2 / 3 ( 1 + M ') and the exponent is greater than unity, 
n oc — £( 1 + 3u ')/ 3 ( 1 + w ). The range of i] corresponding to that of t (0 < t < oo) is 
— oo < 7] < 0, so that the conformal diagram is the lower half of the conformal diagram of 
Minkowski spacetime with the big-bang singularity. Since there is a past null singularity, 
dn is infinite in this well. 

D: w = — 1. Since a = exp(Ht) (— oo < t < oo), n oc — exp(— Ht). Again, the 
conformal diagram is the lower half of the conformal diagram of Minkowski spacetime 
but without the big bang singularity. 

E: w < -1. Since a = (-t) 2 / 3( * 1+w \ rj oc -(-t)( 1 + 3 '«)/ 3 ( 1 +«'). The range of n 
corresponding to that of t (— oo < t < 0) is — oo < n < 0, so that the conformal 
diagram is again the lower half of the conformal diagram of Minkowski spacetime but 
this time with the big rip singularity at the future spacelike infinity. Since the past null 
singularity is null, dn is infinite in this well. 

These results are shown in Fig. |3J For w < —1/3, since the conformal time is 
bounded above, there is a future cosmological event horizon HI] . The light rays 
emitted beyond the horizon never reach an observer at r = 0. Hence, the asymptotic 
region of spacetime cannot be measured, and there is no S-matrix ^l]-§ The proper 
radius of the horizon R c is given by 



for -1 < w < -1/3. R c = H- 1 for w = -1, and R c = (3(1 + w)/(l + 3w))(-t) for 
w < — 1. (Note that the range of t is bounded above for w < —1.) The size of the 

§ However, the situation is not improved much for w > —1/3. see |2U|. 




(12) 




(13) 
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Figure 4. Conformal diagrams of flat cosmological models with dark energy but 
without matter for w > —1/3, w = —1/3, —1 < w < —1/3, w = —1, w < —1. The thin 
solid line is the future event horizon. 

horizon grows for — 1 < w < — 1/3 and remains constant for w — — 1 but decreases for 
w < — 1 because phantom matter violates the dominant energy condition. 

4-2. Including Matter 

When matter is included, lower half parts of the diagrams are replaced with big 
bang singularity. We plot the corresponding diagrams in Fig. Because matter 
was dominated in the past, ] the size of the particle horizon becomes finite even for 
w < —1/3. However, the presence of the cosmological event horizon is not affected by 
including matter, since the universe would be dominated by dark energy in the future 
(see [T3J Uni H3 d] for the detailed discussion of the evolution of the horizon). 
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big bang singularity 



>_1/3 -1<w<-1/3 w<-1 

Figure 5. Conformal diagrams of flat cosmological models with matter and dark 
energy for w > —1/3, — 1 < w < —1/3, w < —1. The thin solid line is the future event 
horizon. 



5. Growth of Structures 

In this section, we study the formation of structures in flat universe models with dark 
energy. First, we consider the evolution of linear matter perturbation, and then on the 
base of it, we consider the nonlinear evolution of perturbation using spherical collapse 
model. 

5.1. Linear perturbation 

The growth of linear matter density perturbation 5 = Spm/pm in various flat 
cosmological models with dark energy is determined by the following equation 



5 + 2HS - AnGp M S = 0. (14) 
Eq. (JT4"j) can be rewritten in terms of g = 5 /a as 

{n M +n x a-^)-p-Al^M + -^-+l(l-w)Q x a'^g = 0.(15) 



dlna? \2 m 2 A J d\na 2 

The exact solution of the linear perturbation equation (JT5j) was obtained by Silveira 
and Waga [221 122] (see also [21] for w — — 1) for a constant w. Denoting the solution as 
D(a)/a, it is given by 

m =2Fl (_±,^i, 1 _± ] jj± a -A t (16) 

a \ 3w 2w 6w Um J 

where %Fi is the Gauss's hypergeometric function. Here D(a) is normalized so that 
D(a) — > a at a — > 0. 

We compute the evolution of density perturbation as a function of scale factor or 
cosmic time. The results are shown in Fig. [HI and Fig. [7| Matter density perturbation 
is normalized by its present value, D = D(a = 1). We assume fl M — 1 — Q x = 0.3. 
The future growth of the density perturbation becomes more suppressed for negatively 
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Figure 6. The linear growth rate as a function of scale factor for w = — 0.3(dot), 
—0.5 (dot-dash), —1 (solid) and — 2(dash) with Qm = 0.3. D(a) is normalized by its 
present value, Do = D(a = 1). 



large w. We find from Eq. ljlfij) that D(a) asymptotically converges toward a constant 
value in the limit of a — > oo for w < —1/3 |25j . 



r(l - £)r(§ - i)r(-i) \n 



i i 



(17) 
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while for > w > — 1/3 



D(a) _ (^.y 3 - a + ^ ) a - , , (18) 



2* 



This is because the time scale for expansion ~ H^ 1 ~ (Gpx) 1 ^ 2 is much less than 
that for growth of density perturbation ~ (Gpm)^ 1 ^ 2 hi the future (px 3> Pm) Elj The 
asymptotic constant in Eq.()17|) and E.f)18[) is the same for both cases and is an increasing 
function of w. Eq.(|17|) and FigElshow that for negatively larger w, the growth of D is 
saturated earlier in the future. This is due to the rapid growth of px relative to pu- The 
difference of D between various w may be suppressed when D is plotted as a function 
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of cosmic time since time elapses more slowly for negatively large w. 



5.2. Spherical collapse 

In the previous section we studied the growth of linear density perturbation. In this 
section we consider spherical collapse model in order to investigate nonlinear growth 
of perturbation. It is already studied both analytically and numerically j2Hj for a 
cosmological constant and the fate of structures for w < — 1 is briefly described in 
[0]. Our approach here will be more analytic than the former and more quantitative 
than the latter for general equation of state of dark energy. 

We consider a evolution of spherical overdense region with uniform density p sp and 
radius R j2Zj- The evolution of radius is described by 

R AtcG AttG ri „ , 

^ = T^-Psp —{l + 3w)px, (19) 

where p sp oc R~ 3 and the mass is M = (4np sp /3)R 3 . The first term represents the 
gravitational force which contracts the sphere, while the second term represents the 
repulsive force for w < —1/3 due to the dark energy and it prevents the contraction. 

Defining dimensionless radius y = R/Ro, where Rq is the radius at present, the 
above equation (fTT)|) is rewritten as, 

(Q M + Q x a-^) -f- 2 ( V -) + i [Q M + (1 - 3w) Q x a~ 3w ] J- (*) 

-I^ + lA Ji M (?)- 2 = 0, (20) 

where A = p sp (Ro) / Pm,o is the ratio of sphere to background density at present and it 
represents nonlinear density contrast. We also define linear density contrast 5q which is 
the density contrast at present if the density perturbation evolve by linear growth rate 
D(a). Then, S is given by 

~p sp (R(a)) 



5q = lim — r-T-Dn = h m — r~: 



- 1 



PM{a) 



D . (21) 



The analytical solution of Eq. (|20j) at a 1 is obtained by 

(22) 



y _ A i/ 3 



lb*™ 



We use the above equation as the boundary condition to solve Eq. (|2T)|) . Either A or So 
can be determined by a condition of y{a = 1) = 1, and hence there is one free parameter. 

In order to relate the mass M of sphere to the density contrast 5q, we assume that 
the overdense region is formed from the la high-density peak of mass fluctuation. Then, 
the linear perturbation at present 5o is equal to the linear mass fluctuation 8m for the 
mass M, 5 = 8m- Here, 5m is given by 

6 M = 2^2 / dkk 2 P(k)W 2 (kr), (23) 
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where P(k) is the power spectrum and W{kr) with r = (2M /VLmHq) 1 ^ is the top-hat 
window function [2E]- We normalise 5m so that 5m = <j 8 = 0.9 at r = 8/i _1 Mpc. We 
assume Qm = 0.3, the baryon density Qb = 0.04 and the Hubble parameter h = 0.7. 

Evaluating the dimensionless radius y(a) in Eq. (|2*U|) for various values of parameters 
Ao, M and 5q, we find that the fate of over-density sphere can be classified into three 
cases. We display these cases in FigJS] The left panel is A — w plane, while the right 
panel is M — w plane (we also show 5 ). 

(i) The region of "Monotonously Expands". For low density region, A < 10, the 
sphere monotonously expands forever. This is because the linear perturbation at 
far future, lim^oo 5(a) = 5 /D lim^oo D(a), is a constant value from Eq. (fT7j) and 
Eq. ()18j) and cannot reach the critical over density 5 C ~ 1.68 (2^] which is 5 when 
the sphere collapses. In the far future, a ^> 1, the radius increases in proportional 
to the scale factor, R oc a. 

(ii) The region of "Collapses" . For high density region, the sphere stops its expansion 
and collapses to R = 0. 

(iii) The region of "Expands- Contracts- ReExpands". For intermediate mass density, 
10 < A < 20 with w < — 1, the sphere stop its expansion and turns around. But 
the repulsive force due to the dark energy, which increases with time, prevents its 
contraction and the sphere re-expands forever. 

In FigEl we show the radius R as a function of time with w = —1.5 for M = 10 14 M Q 
(dot-dash), 4 x 1O 13 M (solid) and 2 x 1O 13 M (dash), as an example of these three 
types. The radius at present is R = [3M/(47rpM,oAo)] 1//3 - As shown in FigEl we n °te 
that the behavior of R is not symmetric about the turn-around time for w ^ —1. 

The fraction of collapsed objects with mass greater than M at the scale factor a is 



where 5m (a) is the mass fluctuation at a, 8m{o) — 8mD(o)/D . The critical density 
5 c (a) is the linear perturbation when the sphere collapses at a. In FiglTUl we show 
F(M, a) as a function of the scale factor with M = 10 n M o (dash), M = 10 13 M o 
(solid), M = 1O 14 M (dot), M = 1O 15 M (dot-dash) and M = 1O 16 M (dot-dot-dash). 
For larger w the objects with the larger mass M > 1O 15 M can form, while for smaller w 
the growth of density fluctuation is suppressed and the mass fraction becomes constant 
in the near future. 

6. Summary 

In this paper, we have classified the future of the universe with dark energy with various 
equation of state. Moreover we have investigated the future structure formation of the 
universe. 



given by |2H] 




(24) 
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Figure 8. The fate of spherical overdense region. Ao is the ratio of spherical density 
to background density at present, M is its mass and do is the linear perturbation at 
present. In the region of "Monotonously Expands", the sphere monotonously expands 
forever. In the region of "Collapses", the sphere stops its expansion and collapses 
to R = 0. In the region of "Expands-Contracts-ReExpands", the sphere stops its 
expansion and turns around. But the repulsive force due to the dark energy, which 
increases with time, prevents its contraction and the sphere re-expands forever. 

We have found that the phase diagrams of the universe in VLm — plane are 
drastically different between w > —1/3 and w < —1/3. Infinitesimally small change of 
w from w = —1/3 to w < —1/3 moves regions of "expands forever" and "recollapses 
eventually" to completely different places. Moreover, for w < —1/3, there is a region of 
"no big bang", which does not exist for w > —1/3. On the other hand, phase diagrams 
of— 1/3 > w > — 1 or w < — 1 are about same, except the "expands forever" region of 
the former is replaced with the "big rip" region of the latter. 

This discontinuity at w = —1/3 can be seen in conformal diagrams too such that 
the event horizon appears for w < —1/3, but w > —1/3. The size of the event horizon 
grows, remains constant and decreases for —1 < w < —1/3, w = —1, and w < — 1, 
respectively. 

Concerning the structure formation, we have found the linear evolution 
asymptotically converges to a constant value for w < and obtained the value in 
the limit of a — > oo as a function of w. 

However, it is not always the case that these linear perturbations eventually 
turn into collapsed objects. We have classified the fate of the overdense regions into 
three categories, i.e., "monotonously expands", "collapses", and "expands-contracts-re 
expands" . And we have shown each region in w — A plane, where A is the ratio of 
density of the spherical region to background density at the present epoch. We have 
also found the largest structure in the future universe as a function of w. 
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Figure 9. The radius R as a function of the cosmic time with w = —1.5 for 
M = 10 14 M Q (dot-dash), 4 x 10 13 M Q (solid) and 2 x 10 13 M Q (dash). 
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